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The status of analysis of processes with hard photon emitted from the initial state (ISR) at 

BABAR is presented. We tag events by the presence of a hard photon in the detector, then 

^ , reconstruct fi'^ ij.~ and several exckusive hadronic final states. The invariant mass of the final 

P^ ' state determines an effective collision center of mass energy at which these measurements 

J^ I can be compared with results from direct e^e~ annihilation process. The first results on 

^"^ . e^e~ — > 2h'^2h~ , where h = tt, K and on e^e~ — » J/ip, obtained with a sample of 89.3 fb~^, 

are presented. 

Measurements of inclusive rj, tv^ , K'^ and p/p production cross sections below the T(45') 
^— s ■ resonance are also presented. These measurements have nearly complete momentum coverage 

""^ — , ' and precision comparable to the best measurements at higher energies. They therefore allow 

^ , for precise tests of QCD predictions and fragmentation models at 10.54 GeV and of their 

scaling properties. 



1 Initial State Radiation processes at BABAR. 



k>( \ Initial state radiation (ISR) processes can be effectively used to me asure e^e annihilation at 

5h I high luminosity e'^'e^ storage rings, such as the B-factory PEP-IP I I I A large mass range 

is accessible in a single experiment, contrary to the case with fixed energy colliders, which are 
optimized only in a limited region. In addition, the broad-band coverage may result also in 
greater control of systematic effects because only one experimental setup is involved. 
The ISR physics program consists mainly on light hadron spectroscopy and measurement of the 
ratio R = a{e^e~ -^ hadrons) / a{e'^ e~ -^ fi^fj.^), which provides the experimental input to 
dispersion integrals for computation of the hadronic contribution to the theoretical estimation of 
the muon magnetic moment anomaly, a^ = [g — 2)^/2 and of the running of the electromagnetic 
coupling to the Z pole, a{M'^). 

The ISR cross section for a particular final state / is related to the cross section for the direct 
annihilation e^e~ — > / through 

^^^ = W{s,x)af{s'), s' = s{l-x)- (1) 

where x = 2ii^*/y^, E* is the radiated photon energy in the nominal center-of-mass frame and 
y^ the nominal cm. energy of the collider. The quantity s' = s(l — x) represents the mass 
squared of the final state system, /. The radiator function VF(s, x) describes the virtual photon 
energy spectrum and can be computed to an accurcy better than 1%. The direction of radiated 
photon is peaked along the initial beams, but for ^/s ~ 10 GeV the fraction at large angle is 
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Figure 1: Left: measured ISR luminosity, integrated over 0.1 GeV for 89.3 fb ^ of BABAR integrated luminosity. 
On the right: the reconstructed j-i'^ fi~ invariant mass distribution in the J/ip region. 



relatively large. It has been shown a 10-15% acceptance for these photons in BABAR. 
The measurement of the corresponding leptonic process e~^e~ — > ^~^^~j provides the ISR lumi- 
nosity. Thus, the Born cross section for a hadronic final state CTf{s') is given by 



AiV^^T, e/ (1 + 5'psj^) 



(2) 



where ^Nj^ {AN^^^^) are the number of detected /7 (^//7) events in the bin of width As' 
centered at s'. ej (e^^) and ^pg^ (^fsr) ^^^ respectively the detection efficiencies and fractions 
of events when the hard photon is emitted by final-state particles. The latter quantity can be 
sizable for the ///i channel, but negligible for most of the low energy hadronic states, which have 
vanishingly small cross sections at the nominal machine energy. 



1.1 The di-muon final state and ISR luminosity at BABAR 

The data used in the analysis presented here correspond to an integrated luminosity of 89.3 fb^ 
collected both at the T{AS) and in the nearby continuum, with the BABAR detector at the 
PEP-II asymmetric i?-factory. 

The BABAR detector is described elsewhere . The information from the tracking system (Sil- 
icon Vertex Tracker and Drift Chamber) is used to measure angles and momenta of charged 
particles. The quartz Cherenkov radiator (DIRC) is the main subsystem for particle identifica- 
tion {KJTx separation is an essential ingredient for these studies). Muons are identified in the 
Resistive Plate Counters installed in the magnet yoke of the BABAR solenoid, while photons 
are detected in the Csl calorimeter. 

Events are tagged by detecting a photon with an energy in the cm. system larger than 3 GeV. 
The invariant mass of the dimuon system determines the effective collision energy. A disadvan- 
tage deriving from the use of ISR events is the invaraiant mass resolution, comparing to the case 
of fixed energy machines. The resolution can be improved by a kinematics fit constraining the 
mass recoiling against the reconstructed final state to be equal to zero, that is the photon mass. 
After the kinematics fit the invariant mass resolution improves from 16 MeV/c^ to 8 MeV/c^ at 
the J/f/' mass. In addition the fit effectively removes background from t^t~ events. Fig. ^left) 
show the derived spectrum of ISR luminosity per 100 MeV energy bins. The BABAR luminosity 
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Figure 2: Left: the four-pion invariant mass distribution. The points indicate the estiimated ISR-type backgrund, 
while the cross-hatched histogram corresponds to the non-ISR background. On the right: e'^e" -^ 27r"'"27r~cross 
section as a function of energy obtained with BABAR ISR data (black dots) in comparison with all previous 

e"'"e~ data. 

of 89 fb~ is equivalent to an e'^e~ machine energy scan in steps of 0.1 GeV, with a luminosity 
integral per point varying from 0.7 pb~ at 1 GeV up to 3.6 pb~ at 4 GeV. This luminosity in- 
tegral provides already a statistically competitive hadron sample, especially in the energy range 
between 1.4 and 3.5 GeV, where very few data are presently available. The systematic error is 
estimated to be of the order of 3% (5% for energies below 1 GeV). 



1.2 Hadronic channels 

A very rich program can be exploited in BABAR from the study of hadronic final states: spec- 
troscopy, form factor measurements, search for exotic states, etc. Besides, measurements of 
exclusive hadronic channels constitute the main approach for measuring R. Currently the ma- 
jor hadronic final states are under study (tttt, KK, 47r, Svr, Gvr, 2K2tt, 4K,pp, KKtt) and more are 
planned. 

Preliminary results are available for the final states with four charged hadrons, namely 27r^27r~, 
2K~^2K~, K~^K~7r~^'jr~ . The discrimination between the three final states is done on the basis 
of the particle identification and on the kinematics fit results for the different mass hypothesis. 
No muon identification was applied since no corresponding background is expected in this 
sample. The background from other hadronic channels is small and is subtracted in each mass 
bin. Overall, more than 70000 events have been selected leading to small statistical uncertain- 
ties. The estimated systematic error is about 5% in the energy region between 1 and 3 GeV, 
dominated by the luminosity determination. Fig.[21^right) shows the derived cross section for the 
4 pions channel in 25 MeV steps. BABAR data are in good agreement with previous available 
results. Moreover, BABAR is the only experiment which cover the entire energy range, with 
an accuracy comparable to the latest precise results from CMD-2'^and SND below 1.4 GeV, 
and much better accuracy than older results from DCI and ADONE above 1.4 GeV. 
The hadronic contribution for this particular channel evaluated using all available e~^e~ data in 
0.56-1.8 GeV energy range is^ al'"^'' = (14.21 ± 0.87exp ± 0.23rad)10-i°, while the r data give 



,hadr 



(12.35 ± 0.96exp ± 0.405(72 )10~^'J. The BABAR data in the same energy region give 



instead^ a^"'^^ = (12.95 ± 0.64exp ± 0.13rad)10 ^°, leading to a substantial improvement. 



1.3 Measurement of J/ip width 

For a narrow state of mass M such as J/V'; decaying into a final state /, the Breit-Wigner can 
be approximated by a 5 function and the production cross section can be written as: 



(j_ 



*J/i/'7^/7 



UTT^TeeBf 



Ms 



W{s,xo) 



(3) 



where xq = 1 — M'^/s. The ratio of /^^/U events from J/ip peak (about 7800 events observed in 
89fb~ , fig. fright)) to the continuum ahows to extract the product!^ 



TiJ/i; 



e e 



B 



j/^-*fi+fi 



(0.330 ± 0.008,iai ± 0.007,y,i)keV. 



(4) 



Background estimation, J/ip hne shape, radiative corrections and Monte Carlo statistics con- 
stitute the main contributions to the systematic error. Using the world average values for the 



branching fraction B 



J/4>~ 



+M+M 



and B 



J/i>- 



(5.61 ± 0.20) keV and Tj/,^ 



^-^-, we derive the J/ip electronic and total widths, 
(94.7 lb 4.4) keV, which represent an improvement 



J/V'— >e e' 

with respect to the present world average values of respectively (5.26±0.37) keV and (87±5) keV. 



2 Inclusive hadron production studies 

BABAR has performed measurements of inclusive production cross sections and fractions of 
etas and charged pions, kaons and protons, both from T(45) decays and in continuum (at 
BcM = 10.54 GeV, below the B-B production threshold). The quark-antiquark colorless system 
created in the continuum events "fragments" into a number of primary hadrons, which then decay 
into stable hadrons. The fragmentation process has been extensively studied at higher energies, 
in particular in the Z^ region and a number of models have been developped and tested. 
Models have been tuned for high energies, but fundamental tests, such as the study of scaling 
properties, need accurate inclusive measurement of hadron productions also at lower energies. 
Previous measurements performed at Argus and Cleo at an energy of ~ 10 GeV were not able 
to cover the full kinematics range. BABAR data are consistent with those results and extend 
coverage up to the highest momentum values with an accuracy at few percent level. 
Fig. 131 shows the differential cross sections in continuum events for 7r,K and p{p) as a function 
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Figure 3: Differential pion (left), kaon (middle) and proton (right) cross section from continuum data in BABAR 
(red points), compared with results by TASSO at cm. energy of 34 GeV (blu triangles) and by SLD at 91.2 GeV 
(green squares). Predictions by JETSET Monte Carlo for different energies are superimposed with corresponding 

color code. 



of Xp, the momentum scaled by half of the center-of-mass energy. BABAR data, obatined from 



a sample of 0.9 fb , are show n a long with those from experiment at higher energies, namely 
TASSoKlat 34 GeV and SLD^! at 91.2 GeV. The nearly fuh kinematics range coverage and 
the comparable precision between BABAR and SLD data, allow to study scaling properties of 
hadronization and to test the predictions by the various fragmentation models. Hadronization 
should be scale invariant. Therefore, cross sections are expected to differ only because of small 
scaling violation effects. At low Xp the difference is due to the mass of the hadrons, with a cut- 
off at Xp ~ 2mh/EcM- At high momentum, substantial scaling violation is expected because 
of running of Og, as it is indeed observed for pions cross-sections. This behavior for pions is 
well reproduced by most of the fragmentation models, as shown for the case of JETSETI^ in 
the same fig. Ol^left). A similar behavior is expected also for the others hadrons, but the data 
show little scaling violation for kaons, and almost no violation for protons. These results, even 
if preliminary, represent already very useful inputs for tuning the simulation of fragmentation 
processes down to an energy of 10 GeV. 



2.1 Test of QCD predictions in the MLLA model 

These data can also be used to test some predictions of QCD in the Modified Leading Logarithm 
Approximation (MLLA/J^. For this purpose, it is convenient to plot the differential cross section 
as a function of the variable ^ = — In(xp). MLLA predicts that a slightly distorted gaussian 
function should be able to fit the data over almost the full kinematic range. Fig. |3 show the 
distributions for vr, K, rj and p/p with the results of the fit. The distorted Gaussian is able to 
describe the BABAR data at the few percent level over the full measured momentum region 
with reasonably small skewness and kurtosis values, consistent with theoretical predictions. 
The position of the peak S,* of the ^ distribution is predicted to decrease exponentially with 
increasing particle mass. The ^* values reported in tabled show a clear disagreement with this 
prediction: the value for pions is quite different from the others, but those for protons and kaons 
are very similar one another, then inconsistent with a continued exponential decrease. This 
qualitative behavior is not peculiar of the BABAR operational energy region, but it has been 
already observed at higher energies and it is particularly evident in measurements performed at 
~ 90 GeV. 



Table 1: Position of the peaks ^* of the ^ distributions for n, k, rj and p{p). 



particle type 


e 


7T+ 


2.36 ± 0.01 


K+ 


1.64 ± 0.01 


T] 


1.44 ± 0.02 


p{p) 


1.61 ± 0.01 



3 Conclusions 



Preliminary results from BABAR on final states produced through ISR demonstrate the high 
physics potential of this sample, which should yield precise measurements of e~^e~ annihiliation 
cross sections. The ratio R = a{e'^e~ -^ hadrons)/ a {e'^e~ -^ fi^fi~), will be measured from 
the sum of exclusive channels, providing input for theoretical determination of the hadronic 
contribution to (^f — 2)^ and a{M'^). The preliminary e^e~ — > 27r"'"27r~ cross section, from 
threshold up to 4.5 GeV has been obtained, with a systematic error of about 5% in the central 
region. The radiative return to J/ijj resonance allows to measure the relative branching fractions 
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Figure 4: DifFerential pion (top), kaon, eta (middle) and proton (bottom) cross section in the continuum data 
(solid symbols) as a function of ^ = — ln(a;p) shown with the result of a Gaussian fit to the data (solid lines) and 
of a distorted Gaussian (G'''^, dashed lines). The distorted gaussian is parametrized as: G'^'''(5; A^, ^, cr, s, fc) = 
— ^=exp \\k + |s(5 — -1(2 + fc)(5 + |s(5 + ^fe(5 j ; where 5 = (^ — ^)/a", ^ is the mean of the ^ distribution, a 
is the square root of its variance, s is the skewness and k the kurtosis. 



with best to date accuracy. 

Preliminary measurements of the inclusive hadronic production cross sections have been per- 
formed at BABAR with precision comparable to the previous measurements at ~ 91 GeV. 
These data will be very valuable for a better comprehension of the fragmentation process and 
in particular to study the scaling properties down to an energy of 10 GeV. The observed discre- 
pancy between data and Monte Carlo predictions on scaling violation for etas, kaons and above 
all protons gives important experimental inputs for tuning models. 
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